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ABSTRACT 


The formulation of a one-dimensional analytical model 
representing steady state operation of a vertically oriented 
gas loaded variable conductance heat pipe is presented. The 
model includes not only the effects of binary mass diffusion 
and axial pipe-wall conduction, but also gravitational 
effects which exist at the vapor-gas interface region when 
the molecular weights of the working fluid and the noncon- 
densible gas differ significantly. Analytical expressions 
for the transport processes (mass, momentum, and energy) 
along with equations of state result in a system of ten 
equations which describe steady state operation of the device. 
Combining these expressions yields a system of two simul- 
taneous highly non-linear differential equations. These 
equations are non-dimensionalized and set up for numerical 
integration using a variable step fourth order Runge Kutta 
method as part of IBM's Continuous System Modeling Program 


(CSMP). 
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cross-sectional area of vapor space 

aXlal cross-sectional area of wick-wall combination 
constant of integration 

control surface 

binary mass diffusion coefficient 

Spatial derivative of Y 

diameter of vapor space 
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Saturation pressure as a function of temperature 
gravitational constant 

a dimensionless function of temperature 
convection heat transfer coefficient 

heat of vaporization 


coefficient of thermal conductivity for wick-wall 
combination 


a dimensionless constant 

length of heat pipe condenser 
condensation mass flux 

molecular weight of noncondensible gas 
pressure 

periphery of outside of pipe 

periphery of vapor space | 

axial conduction heat flux 


power level (heat rejection rate) 





Cle. 


XS 
ede 
eZ 
ae: 


Y4 


maeemwOor axtal Conduceion in wick-wall combination 
universal gas constant 


working fluid saturation density as a function of 
PeMmpecracure 


temperature (absolute) 

velocity 

axial position 

dimensionless length (independent variable) 
dimensionless temperature (dependent variable) 


in We (dependent variable) 


Spatial derivative of Yl (Y3 i) 


Spatial derivative of Y2 (Y4 Dy2) 


volume 


Greek symbols 


u - viscosity 

fe) - density 

T - viscous shear force 
W - mass fraction 
SMDSCEIpES 

a - working fluid 

b - moncondensible gas 
e =e eVaAlOLraceon 

© - ambient conditions 
O - value at x = 0 





fee ROUGH LON 


A. BACKGROUND 

The heat pipe is a heat transfer device capable of trans- 
ferring large quantities of heat under nearly isothermal 
operating conditions. This device is basically a ligquid-vapor 
two phase system which transfers heat by evaporation and 
condensation of a working fluid at opposing ends of a tube. 
There are no moving parts in a heat pipe. To sustain opera- 
tion the condensed working fluid is transferred from the 
condenser to the evaporator section of the device by capillary 
pumping action through a wick structure lining the heat pipe 


inner wall. 


B. CONVENTIONAL HEAT PIPE PRINCIPLES 

The conventional heat pipe in its simplest form consists 
of a hollow tube with wicking material lining its inner sur- 
paeewaowcnhown 1 Figure I. This device is evacuated, then 
filled with an appropriate amount of a working fluid. 

In operation, heat at a given rate from a source 1s ab- 
sorbed at one end of the pipe evaporating the fluid from the 
wick. The vapor then flows along the pipe to the condenser 
end where heat is rejected to a sink and the vapor is condensed 
into the wick. Circulation of the working fluid is maintained 
by a capillary head in the wick which must exactly balance 


pressure drops in the liquid and vapor, as well as any body 
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CONDENSER 
SECTION 


vapor 
wick 


pipe wall 


ADIABATIC SECTION 


EVAPORATOR 
SECTION 


Figure 1. Conventional Heat Pipe 


Pel: 








forces present. The heat load alone controls the satura- 
tion temperature and pressure of the working fluid, which 
become the operating temperature and pressure of the pipe. 
Since the transfer of heat is effected by adding and 
removing the heat of vaporization at saturation conditions, 
the heat pipe is nearly isothermal along its length. How- 
ever, a change in the heat load will automatically change 
the operating temperature (and pressure) of the pipe, given 
a constant sink temperature. This characteristic is unde- 
Sirable for many applications where a constant source tem- 


perature is required. 


C. VARIABLE CONDUCTANCE HEAT PIPE PRINCIPLES 
Ll. General 

If a noncondensible gas is introduced into a heat 
pipe along with the working fluid, the heat pipe becomes a 
passively controlled device which can operate at various heat 
loads while maintaining a more nearly constant source 
memocriatLure . 

During operation the working fluid behaves as described 
above; however, the noncondensible gas is entrained in the 
vapor flow and forms a gas plug at the condenser end of the 
pipe. This plug acts as a diffusion barrier to the flowing 
vapor so that the effective condenser length, and therefore 
condenser heat transfer area, is reduced. For a given heat 
LAG and noncondensible gas loading, the pipe will operate 


at a specific temperature. If the heat load is increased the 


KZ 





saturation temperature and partial pressure of the working 
Fluid, as well as total pressure in the pipe, tend to in- 
crease. In response, the gas plug is compressed which 
increases the effective condenser heat transfer area, thus 
reducing the system temperature rise. A thorough presentation 
on the theory and design of variable conductance heat pipes 
has been provided by Marcus [1]. 
wee lat Pont Theory 

Bienert [2] developed the earliest model used to 
examine the operating characteristics of a gas loaded heat 
pipe. This flat front theory assumed that the noncondensible 
gas was in the form of a plug with a sharply defined vapor- 
gas interface perpendicular to the axis of the pipe. This 
Simple model neglected axial conduction in the heat pipe wall 
and the diffusion of vapor into the stagnant gas plug. MThis 
early model provided much insight into the operating charac- 
teristics of gas loaded heat pipes. However, experimental 
results [1] have shown the model to be inaccurate, primarily 
because it neglects axial conduction and diffusion. 

3. Diffuse Front Theory 

A more accurate model of the gas loaded heat fipe was 
provided by Edwards and Marcus [3]. The model, referred to 
as the 'diffuse front theory', includes effects of binary 
diffusion between the working fluid and the noncondensible gas 
as well as axial conduction in the pipe wall. Marcus reported 
[1] good correlation between experimental results and a com- 


puter program based upon the diffuse front theory. 
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4. Gravitational Effects 

Neither the flat front theory nor the diffuse front 
theory included body force effects on the vapor-gas inter- 
face region. This simplification was probably made because 
much of the early work with heat pipes was related to space 
program applications. 

Kelleher and Batts [4] reported convincing experimen- 
tal results indicating that gravity effects can significantly 
distort the axial temperature profile when the working fluid 
and noncondensible gas are of significantly different molecular 


weights. 


De Cove ecriveE 

The objective of this research was to develop a one-dimen- 
Sional analytical model describing steady state operation of 
a gas loaded variable conductance heat pipe. The model was 
Pomlaeruge not only the effects of axial conduction and binary 
mass diffusion, but also gravitational effects on the vapor- 
gas interface region. An additional purpose of this effort 
was to transform the analytical equations into a set of non- 
dimensionalized simultaneous differential equations which 
could readily be solved numerically using IBM's Continuous 


System Modeling Program III [5]. 
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fT. fORMULATION OF ANALYTIC EQUATIONS 


In formulating the analytic equations describing steady 
state heat pipe operation the differential control volume 
in Figure 2 was developed to represent a section of the heat 


pipe condenser. As implied in Figure 2, only vertically 


(0. = Pp) gAV 
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Pigure Z. Differential Volume Element 
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Oriented heat pipes are considered so that the transport 


processes can be treated as one dimensional. 


A. SIMPLIFYING ASSUMPTIONS 
The following assumptions were made in formulating an 
analytical model using the various transport processes: 

1. The heat, mass, and momentum transport processes were 
assumed to be one-dimensional in the axial direction within 
the vapor space. This assumption restricted this analytical 
model to the specific case of a vertically oriented heat pipe. 
In any other orientation, gravitational effects should be 
expected to cause highly three dimensional spatial variations 
when the molecular weights of the working fluid and noncon- 
densible gas differ substantially. This effect was demon- 
strated experimentally by Kelleher and Batts [4]. 

2. The working fluid and the noncondensible gas were treated 
as thermodynamically independent substances. 

3. Steady state conditions were assumed. 

4. The working fluid was considered to be at saturation 
conditions at all times. 

5. Thermal resistance of the wick-wall structure, in the 
radial direction, was considered to be negligible compared 
to the thermal resistance outside of the pipe. 


6. The noncondensible gas was treated as a perfect gas. 
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pee CONSERVATION OF MASS 
1. Species ‘a' (working fluid) 
Referring to the differential volume element in 


Figure 2 and considering the system on a unit time basis, 
mass of '‘a" in = mass of ‘a' out 


en, 


= qd : 
pvac = Lo V5 + ayl avn) dxlAc +m Pvdx (Gigiee |.) 


which reduces to 
r dq 
Oo = m Pv + Ac =-(0 v_) Gri. 24 


Fick's first law of binary diffusion, as developed 
in reference 6, was used to relate the mass flux of working 
ie oleh Oe to the mass flux of the mixture (pov). For 
the one-dimensional case the scalar form of Fick's law can 


be written, 


BR, Maer = ope. == (ees) 


where the mass fraction of species ‘'a' is defined as, 
ho ee ere 2) 


and the density of the mixture is given by 


Ly 





“ Py Cris: 


Since species ‘'a' and species 'b' occupy the same volume as 


independent substances. 


Then equation II.2 becomes, 


Pv d at 


d . 
aa a! > a> a) = CO (II. 6) 


which represents the conservation of mass for species ‘a'. 


2. Species ‘'b' (noncondensible gas) 


A Similar equation was written for species 'b' noting 


that since species 'b' is a noncondensible gas, Th, SOK 
Then, 
S=(w, ov) = S— (0D, “*b, So (gS 7) 
where the mass fraction of species 'b' is defined as, 
Ww, = Oe, /0 (ieee) 


3. Mass Come enue y; of the Mixture 


Adding the equations for species ‘a' and species 'b' 
(lil@6 and IIlI.7) gave the continuity equation for the 


mixture, 


dw dw 


a d d a b 
ae qx (eviw, +4) ] Pelee wae qa! 0 (Cie 9") 
where D = Dob = Ds is used for simplicity of notation. 
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Recalling equations I1I.4, I1I.5, and II.8, 


which reduced equation II.9 for mass continuity of the mix- 


ture to, 


A Py + S—( pv) a0 (Ir.10) 


C. CONSERVATION OF MOMENTUM 
i ese Ormoe@celes a. aad Species 7b" 
Conservation of linear momentum for the mixture, with 
steady state conditions, was written for the vapor space 


differential element of Figure 3 as, 


ig -* 
= { v(pv + aA) 
cs 
where, dA is the outward directed surface area element and, 
F is the sum of all external forces acting on the fluid in 
the control volume. For the one-dimensional case FP = Be 


The sum of the body forces and surface forces was seen from 


Figure 3 to be, 
ie, = 0g Ac dx - t Fv dx - dp Ac Crie.Ly) 


and the surface integral, representing the net efflux of 


linear momentum in the axial direction from the control 
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Figure 3. Vapor Space Differential Element 


20 





« 
% 


volume, was expanded to give: 


—_ 
tl 


- Sone + (vtdv)(oev+d(pov) )Ac + v m Pvdx 
cs,x aa. 


Acovdv + Acvd(opv) + vm Pvdx (it 2) 


equating equations II.1l1 and II.12, and rearranging yielded, 


Pv dp dv at pv) rn eee 


(O5e eee Ac a anaes = OS we dx ee ao aie 13) 





For steady state conditions the noncondensible gas remains 
stationary with respect to the coordinates fixed to the heat 
Pipe. Therefore the net flux of species 'b' is identically 
zero (1l.e., PVE, = Q since he =0) 5 Noting that 


ov = PV, t+ OMe and Vy, = 50s 


OF 


V 
V = aa WP, = ae 
a a 


Replacing Vea in equation II.13 and rearranging resulted in 


the momentum equation for the mixture, 


m 
d pi Pp P 
eee ie PY ) 55 GS ~ xe —- 
el 


on 





where the viscous shear force t is some function of velocity. 


For example, T = =u v, 1f Poiseuille flow is assumed. 
Zeno weecres ‘Db’ 


Since the velocity of species 'b' is zero, the con- 
servation of momentum equation for species 'b' could be 


written directly by inspection as 


dp, 
ce) Sse (II.15) 
A separate momentum equation for species ‘a’, similar 


in form to equation II.14 could be written. However the addi- 
tional equation would not be useful since it would not be 
linearly independent of the two momentum equations already 


presented. 


D. ENERGY BALANCE FOR PIPE WALL 

In developing a heat pipe energy balance the wick and 
pipe wall were considered together as a composite. Recall 
that thermal resistance of the wick-wall structure, in the 
radial direction, was assumed to be negligible. Figure 4 
then represents an appropriate differential element upon 
which to develop the energy balance. 


By the conservation of energy (referring to Figure 4), 


shia Our 


j= men, 


2iz 





pipe wall 





To 
2 ES eee a eS Cis ree 
x+dx 
Oc + dOc 
Figure 4. Differential Element for Energy Balance 
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Qc cane Neg PY wee Ve rn Ps (l= Tl dx (ii 216) 


but, by Fourier's law of conduction, 


Qc = - kAp (say) 


macorporating equation I1I.17 into equation II.16 and 


rearranging yields the energy balance equation, 


2 h By, 
Geile oes 7 fg : = 
a) k Ap tet) + Ae a 0 (Gis!) 





where, 
(a) the coefficient of thermal conductivity (k) 
must be selected to represent the composite 
Wick-pipe wall structure, and 
(b) the heat transfer coefficient, h, is some 
function of temperature. The correlation 
selected for h(T) depends upon the external 


COO Lind, Come. Ons. 


Ee eOUATIONS OF STATE 
1. Species ‘a' 
Since saturation conditions were assumed to exist 
in the working fluid at all times, tabulated saturation 
properties of the working fluid can be used to generate 


functions for partial pressure, partial density and heat of 
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woeoricaatton of the working Fluid as functions of local 


temperature. Then, 


p = Fa(T) (ae. 2891) 


p = Sa(T) (ito 0) 


The functions Fa(T), Sa(T), as well as Hee repre- 
sent actual tabulated saturation properties of the working 
fluid which can be readily used in a computerized interpo- 
lation scheme. Alternatively Fa(T), Sa(T), and neg!) Gouibe 
be approximated by polynomials or other anlaytic functions 
selected to fit the working fluid saturation curves. 


22 Species 'b' 


As a perfect gas, 


R 
Py = Pp m7 Clie Zak) 
Spee excure., PrOpereLes 
For independent substances, 
and, 
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Although equation II.22 is valid based upon Dalton's law 


of partial pressures, 


only since it 1s not linearly independent of equations 


it is included here for completeness 


II.14 and II.15 developed from momentum considerations. 


Ee SUMMARY OF EQUATIONS 


The following set of ten independent equations in ten 


unknowns describes steady state operation of a vertically 


oriented gas loaded heat pipe. 


ms = + S— (ov) =) 9 
ge = og - Stov4) - Rv a - 
dp, 7 
EBS 2 
ee Pn cron) + £9 BY 
ax 
P, = Py 7 
wi, = e,/P 
p = Pp, t+ Py 
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The ten dependent variables (mo, ie a? Pee Pye 


Ws T) are all functions of axial position. 
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Til. TRANSFORMATION TO COMPUTER MODEL 


A. ELIMINATION OF VARIABLES 

The ultimate goal was to obtain two nondimensionalized 
differential equations with nondimensional temperature and 
mass fraction of species 'a' as the sole dependent variables. 
By reducing the number of equations and dependent variables 
analytically the number of computations required for each 
numerical integration step is significantly reduced. This 
approach was taken to enhance the numerical accuracy and 
Stability of computer results for the model. With this 
objective in mind, the set of ten equations in ten unknowns 
summarizedin Section II.F was reduced to two equations in 
two unknowns by the following steps. 

1. Total pressure (p) appeared in equation I1.14 only and 
therefore could not be eliminated as a dependent variable by 
back substitution. Since total pressure as a function of 
axial position was not of particular interest, p was eliminated 
by assuming that the axial pressure gradient of total pressure 
was negligible (dp/dx = 0). Since the axial pressure gradient 
of total pressure iS primarily created by the hydrostatic 
head of the vapor-gas mixture, this common assumption is 
reasonable. This assumption means that total pressure is 
assumed to remain essentially constant throughout the vapor- 
gas space of the pipe for a given operating power, gas load, 


operating temperature and sink temperature. 
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It 1S Significant to note that although the axial gradi- 
ent of total pressure is negligible, axial variations in 
the partial pressures of both working fluid vapor and noncon- 
densible gas are Significant. For example, using methanol 
as the working fluid at typical operating conditions results 
in an axial partial pressure variation of several atmospheres. 

This assumption eliminated equation II.14 and the variable 
oF 

2. The following simple substitutions were then sequen- 
tially made in the remaining equations. 
a. Equation II.20 was used to eliminate po: 


Sat > Sa(T) 


b. Equation II.19 was used to eliminate Po: 


c. Equation II.21 was used to eliminate Py: 


Bee 
b M. 


Py, T 


ad. Equation II.23 was used to eliminate Pp: 


e. Equation II.4 was used to eliminate op. 


Pa sé Gb (TT) 


9 = Sor ani 


a a 


Zo 





3. The surviving four equations in four unknowns (mh, 


v,T,w_) were: 








du) 

SU Ey d 2G Seb a = 
ale * qx‘54v) dx | ™ a) = ° oan 

Reh i d ,Sa ie 

ne © ax an? at io a 
d i RT, l 
pee a oe = Sa (—— 1)g (IIT. 3) 

a a 

aor | Ps Dip op ) paccnaen = QO (oli sea) 
axe=COK AP 00 k Ap a ; 


4. Equation III.2 was then used to eliminate ms in 


equations III.1 and III.4. 


m = _Acd Sav 
a Py dx We 


) 





Equation III.1 then became: 


dw 
ae! (32 D a =a 6 


Ge Sa” V d d 
ax Wo dx 


ax uw PS ara oc 





which was integrated directly to yield, 


where C is a constant of integration. However, at x = 0 


30 





v = 0 and Sa/w. = 0 


then 


dw 
a 


Cc = ~ PoP “dx |. (ie ros) 





where x = 0 was assumed to be at the condenser-reservoir end 
dw 


__ 
dx 


will be discussed later. 


of the pipe. However, = 0 is a boundary condition as 


Therefore C = 0 


The remaining system of three equations was then, 








du 
a ee oe > a 
Sa v(l — . aaa 0 gs a 6) 
a 
d iL R = 1. 
pe Aol ar = Me i) = Sa (“= Eg ieee) 
a a 
2 fa Ac 
Sly EE Ps h is adayoa Vv _ 
eR Ap’) T=) ~ KAP ax’ u,)? * 9 a 


5. Equation III.6 was finally solved for velocity and used 
to eliminate v in equation III.7 to give a system of two 


equations in T and woe 





a iL R 7 1 _ 
qx Sal = 1)5- 2 = Sa(—— eg erie: 
Jo) a 
2 h Ac du 
dT _ Ps Dip? (es _fg 7d );_DS5a_ a = 0 Cees) 
ax? k Ap oo k Ap dx w (wont) dx 


en 





To solve these two equations simultaneously it was advan- 
tageous to expand and differentiate equation III.3 and to 
expand and rearrange equation III.8 thereby generating a 
system of two second order differential equations in w 
and 2 

After differentiating equation III.3 a considerable amount 


Of algebra led to the following form: 











dw 
ot(= 1) (tT S2+say} - —8(£23) 
we 2 2 
dx a Gicc W 
a 
2 dw 
_ = ve dSa deSanyan. ec 2 at a dS a 
ee cate send) ax) 2 dx dx’ \" “ar °°” 
a ar? W 
a 
| dw 
TSa nel 1 dSa dT Sa Zl 
ma “ay? a ras i) (aa ABS me ee igri) 
a a 


Equation III.8 was manipulated into a similar convenient 


format as: 





2 
a-r ; d tae: he D Sa —_ Ps = , 
dx dx Ap k (wo-l)w. k A 
Ach dw dw 
i a a 1 ik 
mantener or 1 
a a a a 
dD dSa, aT 
+ (Sa aq + D aT Ae (cr. 6) 


Equations III.9 and III.10 represent the final form of the 


desired equations in a dimensional format. 


a2 





B. BOUNDARY CONDITIONS 

To numerically integrate these second order differential 
equations four boundary conditions must be established at 
x = 0. To avoid the obvious singularity when wo, = 0 (at 
the evaporator end of the condenser) integration must proceed 
from the condenser-reservoir end toward the evaporator. 
Therefore x = 0 is chosen to be the condenser-reservoir end 
of the heat pipe. 


Appropriate boundary conditions are then as follows. 


Poe at x = 0 (ig ah 5 thal) 
a, ee Se 0 Eure 12) 
dx 
dw 
$e ike = 0 at x = 0 (ii 1 3) 
4. w = 4) > 0 she We (6. (ia ea 4) 
a oe 


The selection of an initial value of We 1s somewhat arbi- 
trary. However, the initial value of we selected, in conjunc- 
tion with the values of ae and T. used in the model, indirectly 
fix the amount of noncondensible gas and the operating power 
level of the heat pipe being simulated for a particular 


working fluid - noncondensible gas pair. 


C. NONDIMENSIONALIZATION OF EQUATIONS 
To nondimensionalize equations III.9 and III.10 the 


following dimensionless variables were defined. 


3s 





XS 


B gat 


where: 


x/L nondimensional length 
(independent variable) 


<7 
nondimensional temperature 


= T(x) 


prescribed evaporator temperature 


ambient temperature 


ln w logarithm of the mass fraction of 
species ‘a' 


For ease of notation and computer programming the first 


and second derivatives of the two dependent variables were 


represented as follows. 


d(Y1) 





y 3 = axs = bY? 
ee) ee 
y4 = axes = DY2 
Z 
Se Cpyelk) . ~ 2 ae) 
— adxs 2 [a 
axs 
2 
= dtoy2) sagged ON OF 2) 
Die dxXs ie 2 
axs 


Each function of temperature was nondimensionalized with 


respect to the function value at evaporator temperature Us) 
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The following nondimensionalized functions of temperature 


were then developed for direct substitution into equations 


mao eand [11.10. 
D ,; . 
Gl = DIT) nondimensional binary mass 
3 dvEtvicl On) COCtr1I1Clente 
ee Ts) 
Ce = ee nee 5 nondimensional temperature 
e gradient of the binary mass 
diffusion coefficient 
De 
G3 = 5 nondimensional heat of 
fg e vaporization 
- Sa : ; : 
G4 == Salt) nondimensional working 
e fluid density 
aU 
G5 = aS nondimensional temperature 
e gradient of working fluid 
density 
a* (sa) ieee)" 
G6 = a a Ge nondimensional curvature of 
al ae working fluid saturation 
density curve 
G7 = ~— nondimensional convection 
= heat transfer coefficient 


Using the nondimensional variables and functions of 
temperature defined above, equations III.9 and III.10 were 
transformed into the following nondimensionalized system of 
two equations in two dimensionless dependent variables (Yl 
and Y2), and their first and second derivatives with respect 


to axial (dimensionless) length. 
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DY3{(1- Exp(Y2)) ((K1Y1+1)G5 + K1G4)} - DYy4((K1Y1+1)G4) 
= (l1-Exp(Y2)) {Y3(¥3(2K1G5+(K1Y1+1)G6)+K2G5) } 


+ Y4{2Y3((K1Y1+1)G5+K1G4) - G4((K1Y1+1) Y4-K2) } Gutiieskss 


DY3(Exp(¥2)-1) -DY4(K3G1G3G4) = (Exp(Y2)-1)K4G7Y1 


G2Y3 wa Eee } 


+ G1G3¥4(K3(G4(-G= - Texctyay=ay) + G5¥3)) (III.16) 


where Kl through K4 are dimensionless constants: 


La oe 
Keer e= »( T, ) 
2 (2 : ) 
K = 
iS ae Le ae 


(iia) 


He even eter eas ct (1) 


oe: Sap (oa) 


Ps h(T_) L¢ 
Coil 5 ice ie ees 


MUS ocr 


) 


In terms of the nondimensional variables, the boundary 


fouclhenense (equations Till. li through III.14) became: 


ll 
© 


Jig ake = § 20 at XS Gis) 


it 
Oo 


an, V2 = Y2 <0) at XS (Gieliice 9!) 
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siA Woe =) DY = 0 at XS = 0 GEE i 32.0) 
4. Y4 = DY2 = 0 at <o=" 0 (Gh ine Za) 
where Y2 = ln Wes As noted previously with the dimensional 


BOGiiMOfNENIs Downdary condition (equation IIIT.14), the initial 
value of We selected is arbitrary. However the initial value 
of Woe in conjunction with the specified values of ue and. I, 
indirectly specifies the quantity of noncondensible gas and 
operating power level of the heat pipe. Since 0 < w_ < 1.0, 


a 


the corresponding value of Y2 must be negative. 


oy 





IV. DISCUSSION AND RECOMMENDATIONS 


Pee UO LSCUSSION 

The Continuous System Modeling Program (CSMP), developed 
by the IBM Company [5], is a straightforward computer package 
available to numerically solve equations Pees and “lon 
The CSMP package includes a number of optional numerical 
integration methods including a variable step fourth-order 
Runge-Kutta method. 

Although equations III.15 and III.16 are highly nonlinear 
in the dependent variables Yl and Y2, this system of equations 
can be considered linear in the second derivatives of Yl and 


Y2. Therefore these equations can be written in the following 


fOrm: 
(COEF1) (DY3) + (COEF2) (DY4) = Fl CIRVE) 
(COEF3) (DY3) + (COEF4) (DY4) = F2 (Oe, 2) 

where, 
COEFl = {(1-Exp(¥2)) ((KlY1+1)G5+K1G4) } Gly) 
COEF2 = -({((K1LY1+1)G4) (IV.4) 
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Fl = (1-Exp(Y2) ) {yY3(Y¥3(2K1G5+(K1Y1+1)G6)+K26G5) } 


+ Y4{2Y3((K1Y1+1) G5+K1G4) -G4((K1Y1+1)Y4 - K2)} 


Gry. 45) 

COEES == (Bxp(Y2)-1) CIS (65) 

COEF4 = -(K3G1G3G4) (ye) 
F2 = (Exp(Y2)-1)K4 G7 Y1+G1 G3 Y4(K3(G4 

. (G2 ¥3 mee sO 2) + G5Y3)) (IV.8) 





eal (Expt 2) 4) 


Using these definitions, the nondimensionalized equations 
representing steady state heat pipe operation may be numerically 
integrated using the variable step fourth-order Runge-Kutta 


option of the CSMP routine as follows. 


il iL 








aa os (Dpiljaxs, — | sree dXS (IV.9) 
0 0 
a iL 
ee ey 2) dxeme = nee. dXs (IV.10) 
0 0 
if ie a2 
yi o= | (wegjcks = (7 SESE e (IV.11) 
0 axs 


Be, 





1 1 .2 
i =e [ (ajeke Sa) SS omer (IV.12) 
0 0 aXS 
where, 
DMI, eye Cony ey 
DY2 = y4 (IV.14) 
and by Cramer's rule, 
_ (FIL COEF4 - F2 COEF2) 
DY3 = TCOEFICOEFA ~ COEF? COEF) ) (Iv.15) 
oy4 = —(F2 COEF1 - F1 COEF 3) ee 


(COELISZEGOEY a= CGOLE2 CORRS.) 


This model can be used to demonstrate gravitational effects 
SyetnecevapOor—-gas interface region by providing axial profiles 
of temperature and the mass concentration of the working 
fluid. When the molecular weight of the noncondensible gas 
is less than that of the working fluid, the axial profiles 
will show a fairly small interface region between vapor and 
gas. However, when the molecular weight of the noncondensible 
gas is substantially greater than that of the working fluid, 
the profiles will show a distorted (spread out) interface 
region. 

Qualitative results can be obtained from this model by 
uSing additional simplifying assumptions and built-in CSMP 


functions as follows: 
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1. Nominal constant values can be selected to represent 
the convection heat transfer coefficient and the binary mass 


diffusion coefficient. For example, 


h(T) = 1.0 BTU/HR/FT* for free convection 

D(T) = 0.4 PT“ /HR for Krypton-Methanol 
Or 

i OMrmo ce lor eHelum-Methanol 


2. To generate values of working fluid density (Sa) and 


heat of vaporization (h a at any temperature, a CSMP function 


1 
generator can be used to provide quadratic interpolation of 
tabulated saturation properties. Tabulated thermodynamic . 
properties of methanol are available in reference 7. The 
CSMp function "SLOPE" can be used to generate values of 
eioay7at and a’ (Sa) /aT* as Functions of temperature. 

The temperature profile generated by this model for a 
given vertical heat pipe using a particular working fluid- 
noncondensible gas pair, with given input values of Ths 
To, cand w (0), represents heat pipe operation at a certain 
power level and total pressure with a specific load of non- 
condensible gas. The power level, total pressure, and mass 


of noncondensible gas represented could be computed as 


follows. 
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l 
Mass of Gas (ib) = 2p“ f{ 2°(1-w_)axs (IV.17) 
m 4 “v 0 i a 





1 
Power Level (BTU/HR) = PvL f mh, he, oxs (IV.18) 
0 
Sede, ele 
Total Pressure (psia) = De ea (Oj PaO) 77 
a b 
ieny = 9) 


B. RECOMMENDATIONS 
Based upon preliminary work using CSMP, numerical results 
from this model could be improved by the following measures. 

1. Analytical expressions are recommended in 1ieu of the 
special CSMP functions to represent working fluid saturation 
properties as functions of temperature. Expressions for the 
first and second derivatives of density as functions of tem- 
perature could then be found analytically and used directly. 
The resulting expressions should yield more accurate values, 
particularly for the second derivative of density, than are 
achievable through multiple use of the CSMP "SLOPE" function 
with tabulated saturation values. 

2. Since the binary mass diffusion coefficient is a fairly 
strong function of temperature, quantitative results in the 
interface region could be improved considerably by use of 
an empirical expression for D(T) instead of a nominal constant 
value. Reference 8 contains appropriate empirical expressions 


for binary mass diffusion coefficients. 
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3. Quantitative results could also be improved somewhat 
by using an appropriate empirical expression for the convec- 
tion heat transfer coefficient as a function of temperature 
instead of a nominal constant value. Reference 9 provides 
a variety of expressions for h(T) depending upon the convec- 
tive condition being simulated. 

4, With the improvements noted above, this model could be 
used in an iterative loop to obtain preliminary quantitative 
design information concerning vertically oriented terrestrial 
heat pipes when the molecular weights of the operating sub- 
stances are substantially different. 

For example, if the working fluid-noncondensible gas pair 
1s specified, along with the required heat source (evaporator) 
temperature range, ambient conditions, and the maximum safe 
total pressure for structural integrity of the pipe, then 
equation IV.19 could be used to determine the value of w (0) 
POmMiacmucr DOUNdary Condition 219.19. Equation IV.17 could 
then be used to numerically compute the required mass of non- 
condensible gas to load depending upon the heat pipe dimen- 
sions. Equation IV.18 would then yield the maximum heat 
rejection rate obtainable for these conditions. If this heat 
rejection rate is not adequate for the design requirements, 
Plewdesmoanen could Varyeany Gf a number of parameters, Such 
as operating temperature, pipe dimensions, total pressure 
allowed, substances used, etc., until an adequate design is 


found by iteration. 
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